"Moth-eaten effect" driven by Pauli blocking, revealed for Cooper pairs 
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We extend the well-known Cooper's problem beyond one pair and study how this dilute limit 
is connected to the many-pair BCS condensate. We find that, all over from the dilute to the 
dense regime of pairs, Pauli blocking induces the same "moth-eaten effect" as the one existing for 
composite boson excitons. This effect makes the average pair binding energy decrease linearly with 
pair number, bringing it, in the standard BCS configuration, to half the single-pair value. This 
proves that, at odds with popular understanding, the BCS gap is far larger than the broken pair 
energy. The increase comes from Pauli blocking between broken and unbroken pairs. Possible link 
between our result and the BEC-BCS crossover is also discussed. 
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The continuous change from the dilute to the dense 
regime of correlated fermion pairs still is an open prob- 
lem. Although this problem initially arose in the con- 
text of the microscopic theory of superconductivity^^, 
its interest was recently renewed by increasing activ- 
ity in ultra-cold atomic gases. The so-called BEC-BCS 
cross-over between the dilute Bose-Einstein condensate 
of molecules built out of two fermion-like atoms and the 
dense surperfluid state of atom pairs, is a current ma- 
jor questional In the dilute regime, similarities between 
two-atom molecules and excitons should allow their de- 
scription through a composite boson many-body formal- 
ism similar to the one we developed for excitons^. At 
large density, however, excitons suffer a Mott transition 
to an electron-hole plasmgP while Cooper pairs evolve 
toward a BCS superconducting condensate. The physics 
of this BEC-BCS crossover has also been shown to have 
some relevance for Cooper pairs in high-T c cuprates^^. 

In this Letter, we present a conceptually trivial but 
yet unveiled continuity between the Cooper's one-pair 
modeP^ and the BCS superconductivity^. We do it by 
extending the Cooper's problem beyond the single pair 
limit. We start with a "frozen" Fermi sea \Fq) of nonin- 
teracting electrons and we increase the number of elec- 
tron pairs, one by one, within a layer above \Fq) where 
the BCS potential acts. By using this approach, we can 
reach the BCS regime^ continuously starting from the 
single pair limit. 

Although, at the present time, such a pair increase 
seems hard to experimentally achieve, the present anal- 
ysis can at least be seen as a gedanken experiment to 
reveal a possible connection between two famous prob- 
lems in order to more deeply understand the role of the 
Pauli exclusion principle in Cooper-paired states. This 
procedure can also be seen as a simple but well-defined 
toy model to shed some complementary light on the BEC- 
BCS crossover problem since, by changing the number of 
pairs, we do change their overlap. 

The extension of the Cooper's model beyond one pair 



faces a major many-body problem: the exact handling 
of the Pauli exclusion principle between a given num- 
ber of composite particles made of fermion pairs. This 
can be the reason for this extension not to have been 
performed yet. As proposed by Bardeen, Cooper and 
Schrieffer (BCSj^, the smartest way to circumvent this 
difficulty is to turn to the grand canonical ensemble be- 
cause the number of fermion pairs is not fixed anymore. 
This procedure however masks the existing continuity be- 
tween the Cooper's problem and the dense BCS regime. 
This probably is one of the reasons for Schrieffer's claim 4 
that the single-pair picture has little meaning in the dense 
BCS regime. 

We here overcome this quite old many-body difficulty. 
To do so, we start with the equations proposed by 
Richardson^ for the N- Cooper-pair energy in the canon- 
ical ensemble and we manage to solve them analytically 
for an arbitrary number of pairs. This is done by extend- 
ing the method we used to solve Richardson's equations 
for just two pairs^. At the present time, our mathemati- 
cal approach is restricted to the dilute limit on the single 
pair scale. This is why the dense limit is here addressed 
by turning to the grand canonical ensemble and by ex- 
tending the BCS formalism to an arbitrary filling of the 
potential layer. This allows us to show that the solution 
of Richardson's equations we have obtained in the dilute 
limit, remains valid in the dense regime. 

The result we find, proves that the average pair bind- 
ing energy linearly decreases with pair number over the 
whole density range. For the standard BCS configura- 
tion with a potential extending symmetrically on both 
sides of the Fermi level \F) for noninteracting electrons 
- a configuration which just corresponds to fill half the 
potential layer - this gives an average pair binding energy 
reduced to half the single pair value. 

The present work also makes crystal clear how this 
happens. Since Pauli blocking is the only way electrons 
paired by the BCS potential "interact" (see Fig.l), the 
decrease of the average binding energy we find, results 
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FIG. 1: Two Cooper pairs cannot interact through the BCS 
potential given in Eq.(2): this would imply p = k, the 2-free- 
pair state then reducing to zero due to Pauli blocking. 



plainly from the decrease of the number of states avail- 
able for building paired states within the potential layer. 
We can visualize this idea by seeing each added new pair 
as a little moth eating one state, the number of "moth- 
eaten" states increasing linearly with pair number. This 
"moth-eaten" effect, which tends to decrease the effect 
for TV as compared to 1, is actually quite standard in 
the many-body physics of excitons - which also are two- 
fermion states. 

Ground state energy of TV pairs 

As stated above, our work implies handling Pauli 
blocking between a large number of composite bosons. 
This is known to be difficult. However, our knowledge 
on excitons tells us that many important features of 
the composite boson many-body physics are already seen 
when going from 1 to 2 pairs: the effect induced by Pauli 
exclusion principle is already present for two pairs, in 
this way making the understanding for TV pairs far eas- 
ier. This is why, the 2-Cooper pair problem seeming to us 
not out of reach, we seriously looked for the ground state 
energy of 2 pairs, with the idea to extend the procedure 
to 3, 4, TV pairs. 

(i) One-pair energy. The energy of an electron pair 
with opposite spins and zero total momentum, has been 
calculated by Cooper^. It reads E\ = 2s p — e c , where 
ep is the Fermi level of the frozen sea \Fq). In the weak 
coupling limit, the single pair binding energy reduces to 



2n e - 2/poV 



(i) 



po is the density of states taken as constant over the 
potential extension SI. Since the purpose of this Letter is 
to show as simply as possible, the unrevealed consequence 
of Pauli blocking in BCS superconductivity, we accept, 
without questioning it, the "reduced" potential used by 
Bardeen, Cooper, and Schrieffer 
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V is the weak potential amplitude (poV <C 1) while 
Wk = 1 in the energy layer ep < £k < ef + above 
\F ). The main advantage of this reduced potential is 
to be exactly solvable. This allows us to evidence the 
unrevealed physics induced by Pauli blocking between 
Cooper pairs in a sharp way. 

(ii) N-pair eigenstates. Forty five years ago, Richard- 
son has derived^ the exact form for the eigenstates of TV 



pairs. Their energies read as £n = Ri + ■■• + Rn where 
R%, Rn are solution of TV algebraic equations. For 
TV = 2, these equations are 
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plus a similar one with 1 changed into 2 - the equa- 
tions for higher TV's reading as Eq.(3) with all possi- 
ble R differences^. Richardson succeeded to recover the 
BCS resullP^ by solving these equations analytically in 
the infinite-TV limit for a half-filled potential. Today, 
these equations are currently approached numerically for 
small superconducting granules with countable number 
of pairs 17 . However, an analytical solution of these equa- 
tions for arbitrary TV and potential has not been given 
yet. 

(Hi) 2-pair ground state energy. These equations ac- 
tually have a small dimensionless parameter which is the 
inverse of the pair number TV C = p e c above which pairs 
start to overlap - this number increasing linearly with 
sample size. By writing these equations in a dimensional 
form in terms of Zi = (Ri — £\)/e c and by performing 
an expansion in 7 = 1/TV C , we founcP^ that, for a weak 
coupling, the two-pair energy reads, at lowest order in 7 
which turns out to be also an expansion in 1/po 
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Nq = po^l being the number of states in the potential 
layer. 

This result shows that Pauli blocking changes the en- 
ergy of two single pairs {2£\) in two ways: It increases 
the free part by 1/po which just is the Fermi level change 
under a one-electron increase - the extra 2 coming from 
spin. It also decreases the correlated part, one state be- 
ing blocked in the 2-pair configuration. A way to better 
achieve this understanding is to rewrite the single pair 
binding energy e c as 



Eq. (4) then reads 
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Comparison between Eqs.(5) and (6) evidences that the 
correlation energy of two pairs is controlled by the num- 
ber of empty states (Nq — 1) in the potential layer, i.e., 
the number of states available to build the paired config- 
uration. 

(iv) N-pair energy in the dilute regime. It is actually 
possible to solve Richardson's equations along the same 
procedure as an expansion in 7, provided that N/N c stays 
small, a restriction which a priori excludes the dense BCS 
regime, but still corresponds to TV arbitrary large since 
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only N/N c matters. The detailed derivation of this exten- 
sion will be presented in the long version of this Letter. 
Let us here give just a sketch of our procedure. 

Following RefP^, we first rewrite sums appearing in 
the Richardson's equations as 
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where I m = 1 — e~ 2m / p ° v . It can then be shown that, 
when the number of pairs is even, N — 2n, the solution 
for the Zi's at the lowest order in 7 is such that 
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Substitution of Eq. (8) into the Richardson's equations 
leads to n equations for a±, a n which read like 
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We now multiply Eq. (9) by a\ and add to similar equa- 
tions for d2, a n . This leads to 



~ h (al + ... + a 2 n ) + n(n - 1/2) 



(10) 



Next, we turn to the sum of Richardson's equations, as 
given by Eq. (7), with two terms kept, namely 
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(11) 



Using Eqs.(10, 11), as well as the definition of I\ and 1%, 
we can find the sum of Zi at lowest order in 7. From it, 
we get the following expression for the energy of A-pair 
state 
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The same formula for £m can be derived for an odd num- 
ber of pairs, although the form of z^s given by Eq. (8) 
is somewhat more complicated. 

Let us now analyze this result. The first term of £ jv is 
equal to twice the sum £f + ( £ F + p\ )+ —+( £ F + ^^) : 
This just is the energy of A free free pairs added to the 
frozen sea \F ). The fact that we do recover the exact 
normal state energy whatever N, can be a surprise be- 
cause Eq.(12) is a priori derived in the small A/A c limit. 
This led us to think that, most probably, the second term 
of £n also stays valid for A larger than N c . 

(v) Energy in the dense regime. It is first remarkable 
to note that the above result exactly matches the BCS 
condensation energy. Indeed, this condensation energy is 
known to be E BC s = 5P0A 2 with A = 2w c e~ 1/p ° y . As 
2lu c = 51 is the potential extension, Ebcs also reads 



E 



1 



BCS 



-2/p V 



Nnec 
2 2 



(13) 



Nq/2 is the pair number for a potential extending sym- 
metrically on both sides of the Fermi level. The BCS 
result can thus be understood as all up and down spin 
electrons pairs in the potential layer form Cooper pairs, 
their binding energy in this iV-pair configuration being 
half the single-pair energy: This is just Eq.(12) extrap- 
olated to half-filling N = N n /2 for N - 1 ~ N. This 
shows that the "moth-eaten" effect - derived in the di- 
lute limit - seems to stay valid in the dense BCS regime, 
where pairs strongly overlap. 

One important characteristic of the average binding 
energy we find in the dilute limit, is its linear decrease 
with pair number. In order to demonstrate the validity 
of this result in the dense regime, we consider fillings 
different from i.e., a potential extension different 

from fj, — uj c and fi+u> c , the chemical potential fj, being, as 
usual for grand canonical ensemble, afterwards adjusted 
to get the electron number. Textbook BCS formalism^ 
then gives the gap equation as 
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An exact solution exists for fi = ef +f2/2. In the case of 
asymmetrical potential with boundaries still large enough 
to have N S> PoA, we can replace sink -1 by an expo- 
nential. Eq.(14) then gives 



A ~ e- 1/poV 2^/(p-e Fa )(n- p + e Fo ) (15) 

It is possible to show that the condensation energy still 
reads as ^poA 2 , with A now given by Eq.(15). This 
yields Ne c (l — N/Nq), which again agrees with Eq.(12). 

It can be of interest to note that, by inserting Eq.(5) 
into Eq.(12), we can rewrite this condensation energy as 
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since N is the number of occupied states in the potential 
layer while (Nq — N) is the number of empty states. This 
N dependence makes the condensation energy maximum 
when the potential acts symmetrically with respect to 
the Fermi level which precisely is the BCS configuration. 

A last - mathematical - result supporting the validity 
of Eq.(12) at large density, is complete filling. To gain 
condensation energy, empty states feeling the potential 
are required. There is none for complete filling. The only 
possible processes then are electron exchanges. These 
are forbidden within Vbcs- Consequently, condensation 
energy must then reduce to zero. This again agrees with 
Eq.(12) for N = Nq. 

Physical consequences of this A^-pair energy 
(i) Continuity between dilute and dense regimes. The 
above discussion shows that the energy of N Cooper pairs 
given in Eq.(12), although obtained by solving Richard- 
son's equations in the dilute limit, remains valid in the 
dense regime. This supports our understanding, reached 
from the exciton many-body physics, that, due to Pauli 
blocking, the average pair binding energy can only de- 
crease when increasing the pair number, whatever the 
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density. It also reveals a deep connection - missed until to 
now - between the Cooper's picture and the BCS regime, 
in spite of the fact that, as often argued, a strong overlap 
between pairs in the dense regime should destroy any link 
with the Cooper's modeP. This disclosed connection can 
have hidden experimental consequences in superconduc- 
tivity because, as revealed from Eq.(12), paired states do 
have two relevant energy scales: the single pair energy 
e c and the excitation gap A. These two quantities essen- 
tially differ by a factor of 2 in the exponent. This factor 
of 2 however is far from being unimportant because, for 
g-l/poV ver y sm all, it makes the order of magnitude of 
these two quantities quite different. Difference between 
the two factors has already been noted and discussed in 
the literature (see, e.g., p. 169 of RefP}. 

(ii) BEC-BCS cross-over. This connection also offers a 
supplementary route to tackle BEC-BCS cross-over. In- 
deed, in Eagles's and Leggett's approaches, the pair over- 
lap is increased by decreasing the potential V while we 
here increase this overlap by increasing N . These two 
procedures however have some important differences: (i) 
By acting on N, the Pauli exclusion principle blocks more 
and more states while this blocking stays constant when 
one changes V at constant potential extension ft. (ii) 
RefsP^ are based on a BCS wave function ansatz ac- 
cepted as accurate in the dense and dilute regimes but 
more questionable along the crossover^. In contrast, we 
here use the exact wave function obtained by Richard- 
son for the ground state energy of N pairs. In spite of 
these differences, the general conclusion of RefP and the 
present letter stays the same: ground state pairs in the 
dilute and dense regimes are not so much different, a 
conclusion at odds with Schrieffer's clairrP. 

(Hi) Excitation gap. Since the average pair binding en- 
ergy decreases over the whole density range, the reader 
most probably stays with one major question: what con- 
trols the gap in the excitation spectrum of superconduc- 
tors? The answer again is Pauli blocking. When a pair 
is broken, the system not only looses its binding energy, 
but all the remaining unbroken pairs have their average 
binding energy decreased: the two free electrons result- 
ing from the Cooper pair broken by a photon, block two 
pair states (the photon momentum being small but not 
exactly zero). The remaining unbroken pairs feel these 
blocked states when trying to construct their correlated 
state. The latter effect increases with the number of un- 
broken pairs to end in the dense regime, by being far 
larger than the broken pair energy. 

Preliminary results show strong indications that when 
N becomes larger than N c , the threshold energy to break 
a pair achieves the same V and N dependences as A. 
Similar result for the gap change from single-pair t o a 
more cooperative regime was actually found in RefsP^I 
within a variational BCS-like approach, this change going 
along a weak singularity^. 

(iv) Superfluid and virtual pairs. We here deal with 
paired states formed out of all the 2N up and down spin 
electrons added in the energy layer where the potential 



acts. These electrons feel the potential; they are corre- 
lated and form the N pairs we consider in this Letter. 
These pairs are the ones which are "condensed" into the 
same quantum-mechanical state in the BCS wave func- 
tion ansatz. Schrieffer calls thenP "superfluid pairs". 

These "superfluid pairs" have to be contrasted with 
what Schrieffer^l calls " virtual pairs" . The latters corre- 
spond to "electrons excited above the Fermi level" \F) 
of the noninteracting electrons. It is of importance to 
note that the concept of "virtual pairs" is physically rel- 
evant in the dense regime only because the Fermi level 
\F) must not be smeared out too much by interactions 
in order to keep some physical meaning. As a result, 
the understanding of the BCS regime in terms of "vir- 
tual pairs" tends to break in an artificial way a possible 
continuity with the dilute limit. 

These "virtual pairs" are th e ones commonly used to 
give a qualitative understanding^E^l to the BCS conden- 
sation energy, when writing it as a pair number multi- 
plied by a pair energy. Indeed, their number, deduced 
from the width of the BCS distribution change, is of the 
order of — po A. This gives a pair energy of the order 
of A, within an irrelevant factor of 2. From it, it is then 
concluded^ that the "pair energy" must be of the order 
of the gap. It is clear that this conclusion fully relies on 
what is chosen as pair number. By instead taking the 
total number of pairs Nq/2 feeling the potential, as we 
here do - this number being the natural pair number of 
the problem - the same BCS condensation energy gives a 
pair energy exactly equal to e c /2 in agreement with Eq. 
(12). 

We wish to stress that, when compared to the under- 
standing based on "virtual pairs", understanding based 
on" superfluid pairs" provide a natural connection be- 
tween the dilute and dense regimes of pairs. Within these 
"superfluid pairs", the large value of the excitation gap 
is due to many-body effects arising from Pauli blocking 
between broken and unbroken pairs, these many-body 
effects definitely having some physical relevance. 

Conclusion 

We have extended the well-known Cooper's model be- 
yond the one-pair configuration and revealed the simple 
link which exists between this model and BCS supercon- 
ductivity. We show that the average pair binding energy 
linearly decreases with pair number. In agreement with 
our understanding of the exciton many-body physics, the 
Pauli exclusion principle induces a "moth-eaten effect" 
on Cooper pairs, unveiled here for the first time. The 
average pair binding energy in the standard BCS con- 
figuration is shown to only be half the single pair value, 
as a result of their mutual Pauli blocking. This makes 
the excitation gap in the dense regime far larger than the 
broken pair energy. This increase is due to the Pauli ex- 
clusion principle induced by many-body effects between 
broken and unbroken pairs. Our work evidences that su- 
perconductors have a hidden second energy scale - the 
average pair binding energy - which, in the weak cou- 
pling limit, is far smaller than the gap. This result should 
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stimulate new experiments in this very old field. Finally, 
to precisely understand how the isolated pair and BCS 
regimes are connected, can be very valuable in a possi- 
ble approach to the BEC-BCS cross-over within a single 
composite boson many-body formalism^!. 
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